Merging of Dirac points in a two-dimensional crystal 
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We study under which general conditions a pair of Dirac points in the electronic spectrum of 
a two-dimensional crystal may merge into a single one. The merging signals a topological transi- 
tion between a semi-metallic phase and a band insulator. We derive a universal Hamiltonian that 
describes the physical properties of the transition, which is controlled by a single parameter, and 
analyze the Landau-level spectrum in its vicinity. This merging may be observed in the organic salt 
a— (BEDT-TTF)2l3 or in an optical lattice of cold atoms simulating deformed graphene. 
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The recent discovery of graphene has stimulated a 
great interest in the physics of the two-dimensional (2D) 
Dirac equation in condensed matter [1]. The electronic 
dispersion relation e(fc) vanishes at the contact points 
between two bands, the so-called Dirac points D and 
— D (up to an arbitrary reciprocal lattice vector), around 
which the electronic spectrum is linear. Due to the par- 
ticular hexagonal symmetry of graphene, the two Dirac 
points are located at the two inequivalent corners K and 
K' of the first Brillouin zone (BZ). However, that the 
Dirac points are located at high-symmetry points in the 
BZ is not a necessary condition, but a rather special case. 
Indeed, a variation of one of the three nearest-neighbor 
hopping parameters make the Dirac points move away 
from the corners K and K' . If the variation is suffi- 
ciently strong, the two Dirac points may even merge into 
a single one, which possesses a very particular disper- 
sion relation - it is linear in one direction while being 
parabolic in the orthogonal one. The merging of Dirac 
points is accompanied by a topological phase transition 
from a semi- metallic to an insulating phase [2-8] . 

Other physical systems, different from graphene and 
its particular lattice structure, exist where Dirac points 
describe the low-energy properties. Recent papers have 
shown that a similar spectrum may arise in an organic 
conductor, the a— (BEDT-TTF)2l3 salt under pressure 
[9, 10]. Furthermore, it has been shown that it is possible 
to observe massless Dirac fermions with cold atoms in 
optical lattices [3, 11, 12], where the motion of the Dirac 
points may be induced by changing the intensity of the 
laser fields [3]. 

In this Letter, we study in a more general manner the 
motion of Dirac points within a two-band model that 
respects time-reversal and inversion symmetry without 
being restricted to a particular lattice geometry. We 
investigate the general conditions for the merging of 
Dirac points into a single one Dq, under variation of the 
nearest-neighbor hopping parameters. It is shown that 
the merging points may only appear in four special points 
of the BZ, all of which are given by half of a reciprocal 
lattice vector Dq = G/2. Furthermore we derive a sin- 
gle effective Hamiltonian that describes the low-energy 



properties of the system in the vicinity of the topologi- 
cal phase transition which accompanies the Dirac-point 
merging. The effective Hamiltonian allows us to study 
the continuous variation of the Landau-level spectrum 
from oc y/ Bn in the semi-metallic to oc Bin + 1/2) in the 
insulating phase, while passing the merging point with 
an unusual [B(n + 1/2)] 2 / 3 dependence [4]. 

We consider a two-band Hamiltonian for a 2D crystal 
with two atoms per unit cell. Quite generally, neglecting 
for the moment the diagonal terms the effect of which 
is discussed at the end of this Letter, the Hamiltonian 
H(k) reads: 



n( - k) - (/*(fc) '^V) 

The off-diagonal coupling is written as 



-ikR„ 



(1) 



(2) 



where the t m „'s are real, a consequence of time-reversal 
symmetry TL{k) = H*(—k), and R mn = ma\ + 1102 are 
vectors of the underlying Bravais lattice. 

If the energy dispersion e(fe) = ±|/(fe)| possesses Dirac 
points _D, they are necessarily located at zero energy, 
f(D) = 0. From the general expression (2), it is obvious 
that these points D come in by pairs: as a consequence 
of time-reversal symmetry, one has f(k) = /*(— fc), and 
thus if D is solution of f(k) = 0, so is — D. Quite gen- 
erally, the position D can be anywhere in the BZ and 
move upon variation of the band parameters t mn . Writ- 
ing k = ±D + q, the function f(k) is then linearly ex- 
panded around ±D as 

f(±D + q) = -iq ■ (^ t mn R mn cos D ■ Rmn) 

mn 

mn Sin D • R rn n) s (3) 



which has the form q ■ (±Vi — if 2), 
Hamiltonian reads: H.±e> = ±«i • qa* 
of the Pauli matrices a x and a y . 



and the linearized 
+ v-2 ■ q<J y in terms 



Now we consider the situation where, upon variation of 
the band parameters, the two Dirac points may approach 
each other and merge into a single point Dq. This hap- 
pens when D = — D modulo a reciprocal lattice vector 
G = pa\ + <ja*,, where a\ and a*, span the reciprocal lat- 
tice. Therefore, the location of this merging point is sim- 
ply Dq = G/2. There are then four possible incquivalent 
points the coordinates of which are Dq = (pa* + qa*,)/2, 
with (p, q) = (0, 0), (1, 0), (0, 1), and (1, 1). The condition 
/(£>o) = J2 m n(^ 1 Y 3m " t mn = 0, where f3 mn = pm + qn, 
defines a manifold in the space of band parameters. As 
we discuss below, this manifold separates a semi-metallic 
phase with two Dirac cones and a band insulator. 

Notice that in the vicinity of the Dq point, / is purely 
imaginary (v® = 0), since sin(G • R mn /2) = 0. Conse- 
quently, to lowest order, the linearized Hamiltonian re- 
duces to H = q-v%a y , where v% = E m J- 1 )'' n "' t ™n- R m»' 
We choose the local reference system such that v\ = cy 
defines the y-direction. In order to account for the disper- 
sion in the local x-direction, we have to expand /(Do + g) 
to second order in q: 



/(Do- 



-g) = -jg-« 2 -^(-i) /3m "t mn (g- R mn ) 2 • (4) 



Keeping the quadratic term in q x , the new Hamiltonian 

may be written as Tto(q) = ■ 2 ~^< jX + cq y cr y where the 
mass to* is defined by 



1 

m* 



Vf-l^ mn+ 4 R 



2 

mn.x ' 



(5) 
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FIG. 1: Evolution of the spectrum when the quantity A is 
varied and chanyes in siyn at the topological transition (arbi- 
trary units). The low-eneryy spectrum stays linear in the q y 
direction. 

changes its sign at the transition. This parameter A 
therefore drives the transition. In the vicinity of Do, the 
Hamiltonian becomes Tt(q) = Ho(q)+Aa x , or explicitly : 



H(q) 



2m* 



icq y 



A 



2m ' 



icq y 



(8) 



where R mn ,x is the component of R mn along the local 
x-axis (perpendicular to v®). The terms of order q^ and 
q x q y are neglected at low energy. The diagonalization of 
Tio(q) is straightforward and the energy spectrum e = 

(q x /2m*) 2 + c 2 q y has a new structure: it is linear in 
one direction and quadratic in the other. From the linear- 
quadratic spectrum which defines a velocity c and a mass 
to*, one may identify a characteristic energy : 



E mn (-i) /3 " 



1 ^mn-^rnri] 



E mn (-i) & » +1 ^ 



H 2 



(G) 



Up to now, we have discussed the merging of the two 
Dirac points from a "dynamical" point of view, following 
their motion in the BZ when varying the band parameters 
until Dq is reached. We now consider the low energy 
Hamiltonian around Do even before the two Dirac points 
coincide. In the neighborhood of the transition when 
/(Do) = 0, there is a finite gap 2| A| at the Do point (see 
Fig. 1), where the quantity 



A 



i) 



(7) 



with the spectrum e = ±w(A + q 2 /2m*) 



c 2 ql 



The 



Hamiltonian (8) has a remarkable structure and describes 
properly the vicinity of the topological transition, as 
shown on Fig. 1. When to* A is negative (we choose 
m* > without loss of generality), the spectrum ex- 
hibits the two Dirac cones and a saddle point in Do (at 
half distance between the two Dirac points). Increas- 
ing A from negative to positive values, the saddle point 
evolves into the hybrid point at the transition (A = 0) 
before a gap 2A > opens. Due to the linear spectrum 
near the Dirac points, the density of states in the semi- 
metallic phase varies as |e| at low energy and exhibits a 
logarithmic divergence ln(||e| — |A||) due to the saddle 
point. At the transition, it varies as yfef and then a gap 
opens for A > [4]. 

The topological character of the transition is displayed 
by the cancelation of the Berry phase at the merging of 
the two Dirac points. The spinorial structure of the wave 
function leads to a Berry phase \ § X78 q ■ dq, where q = 

arctan - b 1 ^/ 1 } . Near each Dirac point, a = arctan— , 
ixef(q) y 9* 

whereas 9 q = arctan 2m q 2 9y near the hybrid D point at 

the transition. Therefore, the Berry phases ±7r around 

each Dirac point annihilate when they merge into Dq [4] . 
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We now turn to the evolution of the spectrum in a 
perpendicular magnetic field B. After the substitution 
Qx ~~ * Qx — e By in the appropriate gauge and the intro- 
duction of the dimensionless gap <5 = A/(m*c 2 u; 2 /2) 1 / 3 , 
in terms of the cyclotron frequency lo c = eB/m*, the 
eigenvalues are e n = ±(A/ S) E n (5) where the E„ are 
solutions of the effective Schrodinger equation 

E n (5)i> = [P 2 + (6 + Y 2 ) 2 - 2Y] 4, = 7Uffi> (9) 

with [Y,P] = i. The effective Hamiltonian is therefore 
of Schrodinger type with a double-well potential when 
S < 0, which becomes the quartic potential Y 4 — 2Y at 
the transition and then acquires a gap for S > 0, with a 
parabolic dispersion at low energy (see Fig. 2). For large 
negative 5, one recovers two independent parabolic wells 
with an energy shift ±2y/S equal to half the cyclotron 
energy. Therefore, as seen in Fig. 2(a), the lowest level 
has zero energy, and the first levels are degenerate : one 
recovers the physics of independent Dirac cones in a mag- 
netic field, and the effective energy levels are given by 
E n = 4nVlJ|. 

V{Y) E„ V(T) E„ 




FIG. 2: Potential profile V(Y) = (S + Y 2 ) 2 -2Y and effective 
energy levels E n for 8 — —4 (a) and S = (b). 

The complete Landau levels spectrum e n (S) is shown 
in Fig. 3. The value of the Hall integer is indicated 
in the gaps between Landau levels. For negative 5, one 
recovers the spectrum of the Dirac cones, with odd values 
of the Hall integers - the absence of even values reflects 
the twofold valley degeneracy of the Landau levels and 
the presence of a zero-energy Landau level. When — 5 
vanishes, approaching the transition, the level degeneracy 
is lifted, and gaps with even Hall integer open. A simple 
WKB analysis of the lowest level shows that it splits as 
to oc ie - ^ 5 ! 37 ~ ±e~#' A 1 l B . The energy levels scale 
as e„ oc B 2 / 3 y/E n (A/B 2 / 3 ), with the following limits 

m* A < , semi-metal — > e n oc ±V nB (10) 
A = , transition -> e„ oc ±[(n + 1/2)B] 2/3 
771* A > , insulator -> e„ = ±[A + #(n + 1/2)B] 

Note the shift n — > n + 1/2, a consequence of the ±7r 
Berry phases. 




FIG. 3: Energy levels e n (5)/(m*c 2 uj 2 /2) 1/3 as a function of 
the dimensionless parameter S oc A/B 2 ^ 3 . The dots on the 
5 = axis indicate the semi-classical levels of the quartic 
Hamiltonian [4]. 

We now consider two specific situations in which the 
merging of Dirac points may be observed. The first exam- 
ple is a variation of the standard graphene tight-binding 
model, where the three hopping integrals between nearest 
carbon atoms are assumed to be different: 

/(As) = too + t w e- lk a ' + t 01 e- lk a - . (11) 
A merging at Do = (p a l + 9 a 2)/2 is possible if 

too + M) p <io + (-l)%i =0 . (12) 

Choosing i o = t' > and t w = t i = t > 0, 
Eq. (12) has a solution (£' = 2t) for p = q = 1, at 
Do = (a* + o,*,)/2, that is at the M point located at the 
edge center of the BZ [4]. Even if the hopping integrals 
may be modified in graphene under uniaxial stress [5, 6], 
it seems impossible to reach physically the merging con- 
dition. An alternative for the observation of Dirac points 
has been proposed with cold atoms in a honeycomb op- 
tical lattice. The latter can be realized with laser beams, 
and by changing the amplitude of the beams, it is possi- 
ble to vary the band parameters and to reach a situation 
where the Dirac points merge [3]. 

The organic conductor a-(BEDT-TTF) 2 I 3 is also a 
good candidate for the observation of merging Dirac 
points. In order to study the low-energy spectrum 
(close to half-filling) , the original description with four 
molecules per unit cell can be reduced to a two-band 
model in a tetragonal lattice, with the following disper- 
sion relation [9, 10, 13]: 

f(k) = ho + t 1Q e ik ai + toie lk - a * + t n e ik -^ +a ^ . (13) 
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In this case, the generic spectrum exhibits two Dirac 
cones the positions of which are given by [5] 




-3-2-10 1 2 3-2 -1 1 2 




-3-2-10123 -3-2-10123 



FIG. 4: Top figures: left) motion of the two Dirac points D 
and —D for the case too = toi = 1, iio = 2 + tn under varia- 
tion oftu; right) phase diagram. The V and X lines separate 
the semi-metallic phase from the band insulator (grey). The 
Dirac points move from X = (0, 1) to the F = (0, 0) point 
when tn varies from to — 2. Bottom figures show the iso- 
energy curves for tu = —1.5 (left) and — 2 (right). 

Upon variation of the band parameters, the two Dirac 
points may merge when 

too + (-l) p tio + (-l)'toi + (-l) p+9 tn = . (15) 

Katayama et al. have considered the situation (in our 
notations) where too = ^oi [9] an d shown the possibility 
of a transition from a massless "Dirac" phase to a gapped 
phase at a hydrostatic pressure ~ 40 kbar [10]. In Fig. 
4, we show the evolution of the Dirac points in the BZ 
(as in Rcf. [9]), and more important, the evolution of the 
spectrum for a particular variation of the band parame- 
ters. The two Dirac points merge at the BZ points T and 
X for special values of the band parameters. 

The scenario can be even richer. One may imagine 
a situation where, when varying a band parameter, the 
Dirac points disappear and then reappear at a different 
Do point of the BZ. In Fig. 5, the Dirac points move 
from X\ to r, where a gap opens. For further variation 



of the band parameter, the gap persists until a new pair 
of Dirac points appears at a different position Xi in the 
BZ, and disappears again at the fourth special point X 3 . 




FIG. 5: Motion of the Dirac points for model (13), with 
too = — 1, tm = —0.5, toi = 1.4, while varying —3 < tu < 3 
(vertical axis). 



We finally consider the effect of non-zero diagonal 
terms in the Hamiltonian (1). When there is inversion 
symmetry, one has H.22(k) = Wu(-fe). Moreover, time- 
reversal symmetry implies that these diagonal matrix el- 
ements are symmetric functions in k, and that their ex- 
pansion near the hybrid point Do = G/2, has no linear 
term. Therefore all considerations discussed above re- 
main valid, although the Dirac and hybrid points are no 
longer necessarily at zero energy. 

In conclusion, we have studied under which general 
conditions the merging of Dirac points may occur, mark- 
ing the transition between a semi-metal and a band in- 
sulator. We have fully described the vicinity of the tran- 
sition by means of an effective 2x2 Hamiltonian. Al- 
though it has been constructed to describe the low en- 
ergy spectrum near Do, this Hamiltonian is appropriate 
to describe both valleys around the D and — D points 
avoiding the use of a 4 x 4 effective Hamiltonian as it is 
usually done. It may even provide an effective description 
of graphene, which could be useful, e.g., in accounting for 
intervalley scattering in a disordered system. 

We recently learned of a related work which proposed 
the existence of hybrid points in the absence of time re- 
versal symmetry in VO2 /T1O2 heterostructures [14] . We 
also became aware of a recent independent work on the 
anisotropic honeycomb lattice in a magnetic field, which 
has some overlap with ours [15]. Finally we wish to men- 
tion a recent paper on mimicking graphene physics with 
ultracold fcrmions in an optical lattice [16]. 
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